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Continuous Quantum Measurement and the Quantum to Classical Transition 
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While ultimately they are described by quantum mechanics, macroscopic mechanical systems are 
nevertheless observed to follow the trajectòries predicted by classical mechanics. Hence, in the regime 
defining macroscopic physics, the trajectòries of the correct classical motion must emerge from 
quantum mechanics, a process referred to as the quantum to classical transition. Extending previous 
work [Bhattacharya, Habib, and Jacobs, Phys. Rev. Lett. 85, 4852 (2000)], here we elucidate this 
transition in some detail, showing that once the measurement processes which affect all macroscopic 
systems are taken into account, quantum mechanics indeed predicts the emergence of classical 
motion. We derive inequalities that describe the parameter regime in which classical motion is 
obtained, and provide numerical examples. We also demonstrate two further important properties 
of the classical limit. First, that múltiple observers all agree on the motion of an object, and second, 
that classical statistical inference may be used to correctly track the classical motion. 

PACS numbers: 03.65.Bz,05.45.Ac,05.45.Pq 



I. INTRODUCTION 

Macroscopic mechanical systems are observed to obey 
classical mechanics to within experimental error. How- 
ever, the àtoms which ultimately make up these systems 
certainly obey quantum mechanics. Therefore, the ques- 
tion of how the observed classical mechanics emerges 
from the underlying quantum mechanics arises immedi- 
ately. This emergence, referred to as the quantum to 
classical transition, is particularly curious in the light of 
the fact that the equations of motion for the trajectòries 
of classical mechanics are nonlinear, and can therefore ex- 
hibit chaos, whereas even a proper quantification of chaos 
in quantum mechanics has been difficult to obtain Q. 

Note that the task of explaining the quantum to clas- 
sical transition (the QCT) is essentially a practical qües- 
tioni it is a question of explaining why real systems, such 
as nonlinear pendulums, baseballs, and other systems 
which can be be built and observed in the laboratory 
obey classical mechanics (at least to within any experi- 
mental error). It is not a question of obtaining classical 
mechanics precisely as a formal limit of quantum me- 
chanics. In fact, due to the absence of chaos in closed 
quantum systems fl , and the non-commutativity of the 
twin limits h — > (the semi-classical limit) and t — > oo 
(the long-time limit necessary to describe chaos) efforts 
to extract classical chaos as a formal limit have been less 
than successful. 

If one describes macroscopic objects sufficiently realis- 
tically using quantum mechanics, then it should be possi- 
ble to predict the (often chaotic) trajectòries of classical 
dynamics. In order to do this, it is important to realize 
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that all real classical systems are subject to interaction 
with their environment. This interaction does at least 
two things. First, it subjects the system to noise and 
damping [j| (as a consequence all real classical sys- 
tems are subject to noise and damping — even if small), 
and second, the environment provides a means by which 
information about the system can be extracted (effec- 
tively continuously if desired) , providing a measurement 
of the system |,|, |,|7H. 

Two levels of description have been used to discuss the 
QCT. The first utilizes the decoherence resulting from 
tracing over the environment to suppress quantum inter- 
ference Q : it is assumed that no dynamical information 
about the individual system has been extracted from its 
environment. In many circumstances this alone can lead 
to an effectively classical evolution of a phase space dis- 
tribution function |Ï0[ . As mentioned above, a more fine- 
grained description is achieved when the environment is 
taken to be a meter that is continuously monitored, lead- 
ing to a 'quantum trajectory unraveling' of the system 
density operator conditioned on the measurement record. 
If one averages over all possible measurement results, the 
description reverts to that at the level of phase space dis- 
tributions. However, the fine-grained description which 
explicitly incorporates monitoring of the environment is 
required to understand the QCT at the level of extracting 
classical trajectòries from the quantum substrate. 

An example of an environment that naturally provides 
a measurement is that of the electromagnètic field which 
surrounds the system. Monitoring this environment con- 
sists of focusing the light which is refiected from the sys- 
tem, allowing the motion to be observed. If the envi- 
ronment is not being monitored, then the evolution is 
simply given by averaging over all the possible motions 
of the system. Classically this means an average over 
any uncertainty in the initial conditions, and over the 
noise realizations. However, in the absence of explicit 
observation (monitoring the environment and recording 
the evolution) it is impossible to obtain classical trajectò- 
ries: the system must be described by an ever broadening 



2 



probability (or pseudo-probability) distribution in phase 
space. This is an experimental truism, and therefore ap- 
plies regardless of whether the system is being treated by 
a classical or quantum mechanical theory. 

Since all classical systems are subject to environmen- 
tal interactions, and since measurement is necessary to 
deduce the trajectòries of classical motion, it may be ex- 
pectcd that such environmental interaction, and the as- 
sociated measurement process, will need to be included 
in a treatment that is adequate enough to predict the 
emergence of classical motion from quantum mechanics. 
Indeed, recent work by a number of authors has made 
it increasingly clear that this provides a natural explana- 
tion for the emergence of classical motion, and, therefore, 
a resolution of thcproblcm of the emergence of classi- 
cal chaos @ H H |H, H- Fortunately, the quantum 
theory of environments and continuous measurement is 
now sufficiently well developed that their effects can be 
treated in a fairly straightforward manner, the emergence 
of classical dynamics verified, and the mechanism of the 
quantum to classical transition elucidated. 

Detailed studies of the QCT are particularly timely be- 
cause current experiments in quantum and atòmic òptics 
and condensed matter physics arc bcginning to probe this 
transition directly in both ensemble and individual sys- 
tem cases [|Ï6| . Our approach here is to present a general 
formalism for understanding the transition: more focused 
analyses appropriate to specific experimental situations 
can easily be developed based on this general approach. 

In the following we examine the QCT in some detail. In 
Section |0], we examine how macroscopic systems may be 
treated, including environmental interactions and mea- 
surement. In Section III, we derive inequalities that de- 



scribe the regime under which classical motion emerges. 
In Section IV, we show that, in addition, classical state- 
estimation will work in the classical limit. In Section [v|, 
we provide two specific numerical examples showing that 
classical motion is indeed obtained in the regime pre- 
dicted in Section |II]|. We finish with some concluding 



remarks in Section VI 



II. 



DESCRIBING THE MOTION OF 
MACROSCOPIC OBJECTS 



A macroscopic object is composed of a very large num- 
ber of quantum degrees of frcedom. For example, we 
can consider the motions of the àtoms which comprise a 
massive object, and these are all coupled together by the 
inter-atomic forces. The equations of classical mechanics 
are supposed to describe the dynamics of macroscopic 
quantities, such as the center of mass; classical motion 
is not observed in the entire many particle phase space. 
Hence, one should consider a change of variables, so as 
to write a Hamiltonian in terms of the center-of-mass 
coordinate, X (with conjugate momentum P). This co- 
ordinate is coupled to all the other coordinates x\ (in 
a sòlid we might refer to these as the internal phonon 



modes, for example). Under the assumption that none of 
these environmental modes is strongly perturbed by the 
dynamics, it is sufficicnt to treat them as harmònic oscil- 
lators, and to couplc them to the center-of-mass motion 
via the Hamiltonian [0, [| 



H = H cm (X,P) + J2 



1 2 i K i i 



X) 



, (1) 



where pi are the momenta of the internal modes, and Kí 
gives the strength of the coupling between the center of 
mass and the internal degrees of freedom. The center of 
mass now constitutes effectively an open system interact- 
ing with an environment consisting of a large number of 
harmònic oscillators with a correspondingly large range 
of freqüències. This is the starting point for a treatment 
of quantum Brownian motion, as developed by Caldeira 
and Leggett Under general conditions, the phonon 
environment can be treated as a heat bath, and in the 
double limit of weak coupling to this bath and high tem- 
perature, it is possible to write a very simple master equa- 
tion for the center-of-mass motion |H|: 



-[H cm (X, P),p]-k env [X, [X,p]} 



(2) 



where k env is determined by the Kí and the temperature. 
This provides a simple description of the effects of the in- 
ternal degrees of freedom upon the macroscopic motion 
of an object for which frictional effects are negligible, and 
the heating due to the noise is not significant over time 
seales of interest. If one wanted to treat damped clas- 
sical systems, then one would relax the weak coupling 
approximation so as to give a master equation that ex- 
plicitly contains damping. However, for simplicity, we 
will restrict our attention here to classical Hamiltonian 
systems. 

Another important environment which we need to con- 
sider is the quantum electromagnètic field. This interacts 
with the object, and provides a natural mechanism for 
measurement of the center-of-mass position X . In gen- 
eral, macroscopic objeets are bathed in light from all di- 
rections, and the light that is reflected may be monitored 
by a large number of observers. Since we are considering 
a one-dimensional system, and since we wish to use the 
simplest description which captures the essential aspeets 
of the measurement process, we restrict ourselves to in- 
teraction with an electromagnètic field in one dimension. 
In particular, we consider a làser reflected from the ob- 
ject such that the phase shift provides information about 
X. Pcrforming an analysis of such a measurement, one 
finds that the evolution of the system, conditioned upon 
the measurement record, may be written as a stochastic 
master equation in the Itò formalism JTÏJ as Jï^, [ïíj 

dp = -{^[H cm {X,P),p] + k([X^X,p} + -2XpX^}dt 
n 



Vk 



{{Xp+ pXÏ) - pTr p(X + X*)}dW 
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= -{j-[H cm (X,P),p} + k[X,[X,p}}}dt 
\fk 

+—{[X,p}+-2pTrpX}dW, (3) 

where the observed measurement record is given by 

dy = Tr pXdt + -=dW. (4) 

In these equations, dW is a white noise generating a 
Wiener process, k gives the strength of the interaction 
between the light and the object and is proportional to 
the power of the làser, whereas k gives the rate at which 
information about the system is obtained. When no in- 
formation is obtained (i.e., k = 0), or if the measurement 
record is averaged over, the stochastic master equation 
(||) reduces to the ordinary master equation (j^). 

The ratio rj = k/8k, called the efficiency of the mea- 
surement (§) , is a measure of the fraction of the refiected 
light that is actually detected by the observer in making 
the measurement. As will be clear from our discussion in 



Sec. III, though both k and k arise from the interaction 
of the system with the measurement environment, they 
play very different roles: whereas h 2 k represents a noise 
on the system that leads to spreading out in phase space, 
k provides information about the system leading to lo- 
calization around individual trajectòries. The fact that 
Eq. (^) leads to a (completely) positive evolution for all 
initial conditions if and only if r) < 1 |^0| is a particu- 
lar case of the general information-disturbance principies 
in quantum mechanics: any process that leads to infor- 
mation about a system must produce at least a minimal 
unavoidablc disturbance. 

If there exist múltiple observers dividing the avail- 
able refiected light up among them, then each sees an 
evolution with a valuc of r\i < 1 (and, for positivity, 
Si Vi ^ Xh w ith a different noise realization for each 
observer EM . This is certainly the case in reality, where 
each observer usually captures only a small fraction of 
the available light. In the regime in which classical mo- 
tion is obtained (which we will refer to as the classical 
limit), all observers must agree on the motion of the sys- 
tem to within experimental error, and we consider this 
question at the end of the next section, and in our nu- 
merical examples. 

Since the form of the equation resulting from interac- 
tion with the infernal modes is the same as that which 
results from failing to monitor the light which is being 
used to probe the system, we can take this environment 
into account in the same way that we take múltiple ob- 
servers into account, that is, by taking an appropriatc 
value of rj < 1. (The measurement constant k is then 
adjusted to include the contribution from k env ). 

The stochastic master equation (||) constitutes our 
description of the evolution of the center-of-mass of a 
macroscopic object. In the following we will show that 
this description, while very simple, is sufnciently realis- 
tic to obtain the correct classical motion in the classical 



limit. It should also be noted that while we have chosen 
to measure the position X, the analysis which follows 
suggests that the extraction of the classical limit is not 
sensitive to the precise observables which are measured; 
as long as the measurement provides sufficient informa- 
tion about the location of the system in phase space, the 
classical limit will be obtained. For example, a contin- 
uous measurement of momentum will suffice, as long as 
the forces on the system are spatially dependent. In fact, 
other authors have provided numerical support for this 
view by showing that quantum state diffusion (using a 
measurement interaction which includes damping) ÍVm 
or a simultaneous measurement of position and momen- 
tum jÏ5| are sufficient to induce the QCT in the same 
manner. 

If the correct classical mechanics is to be obtained, 
two conditions need to be satisfied. First, it must be 
possible to observe the system so that its center of mass 
(and all other degrees of freedom considered classical) 
is known sufnciently accurately on the scale of the po- 
tential and relevant dynamical timescales. Second, these 
observed vàlues, which we might identify as noisy coun- 
terparts of the means of the sufnciently well localized 
distribution, x = (X) and p = (P), should evolve ac- 
cording to the classical Hamiltonian H cnï (x,p) which has 
the same functional form as the quantum Hamiltonian 
ií cm (I, P), with deviations small compared to the clas- 
sical seales. In other words, the existence of the quantum 
to classical transition implies that in the classical limit 
we can replace the quantum operators with effective clas- 
sical dynamical variables, 



III. 



H(X,P)^H(x,p). 



INEQUALITIES GOVERNIN G THE 
CLASSICAL LIMIT 



(5) 



We now ask for the parameter regime in which the evo- 
lution reduces to classical motion. As explained in the 
previous section this means that the quantum distribu- 
tion remains sufnciently localized (such that the system 
can be said to be executing a trajectory), and that this 
trajectory, characterized by x = (X) andp = (P), follows 
that of the classical motion, generated by H cm (x,p). 

We proceed by first writing down the equations of mo- 
tion for the first and second moments of X and P. From 
Eq. (0) these become 



dx = —dt+VkV x dW, 
m 

dp = (F(X))dt+V%C xp dW , 



(6) 
(7) 



and 



dV x = 



-C XTI — kV 2 



dt + VhK xxx dW , (8) 



dV p = [2h 2 k - kC 2 xp + 2d x FC xp ] dt 
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+d 2 x FK xxp dt + V%K xxp dW , 



dCxp — 



-V p - kV x C xp + 8 X FV 3 



d.t 



-\d 2 x FK xxx dt + \[ÏK xpp dW , 



where 



V x = (X 2 ) (X) , 



v p = 



(P z ) - (P) 



Cxp = -(XP + PX)-(X)(P) 



(9) 



(10) 



(11) 

(12) 
(13) 



are the second cumulants, and the K's are the third cu- 
mulants defined by 



K. 



bc 



ABC 



(: AB :){C) 
:AC :){B)- 



-(A)(: BC :> 
2(A)(B)(C), 



where A, B, C can be X or P, and the colons denote Weyl 
ordering of the operator products. In the above equa- 
tions we use the simplified notation F = F(x) = F((X)) 
and expand F in a Taylor series about X = x truncated 
to second order. Without this truncation higher deriva- 
tives of F would appear in the equations for V p and C xp , 
multiplied by higher powers of the widths or by higher 
cumulants. Truncating the power series for F in this way 
is a good approximation so long as the distribution is suf- 
ficiently localized about x and p. Examining Eq. (^]), one 
sees that to maintain classical motion for x one needs 
(F(X)) « F(x), which happens when the system is lo- 
calized enough so that V x d x F(x) <C 2F(x). It is the task 
of the measurement to maintain such localization, and 
numerical studies [Q have shown that it can indeed do 
so. 

At this point, it is perhaps instructive to look at the 
origin of localization of the individual trajectòries. The 
density matrix obtained by solving Eq. (^|) is conditioned 
on the measurement record Eq. (||), or equivalently, by 
the noise realization dW. Averaging over these realiza- 
tions results in the density matrix of the unobserved sys- 
tem which can also be obtained by solving Eq. The 
second cumulants cr 2 x , cr xp , and a px of that distribution 
are related to the corresponding cumulants for each tra- 
jectory by the relations: 



= (V x )w + yaxw(x,x) , 
= (C xp ) w + var w (x,p) . 



&PP = (V P )w + vai w (p,p), 



Xp 



(14) 



where (-)w an d varvy(-,-) represent respectively the 
means and (co-)variances of the quantities, when con- 
sidered as distributions over trajectòries. The Wiener 
process damps the first term on the right hand side of 
each equation by a term proportional to —k, at the same 
time compensating this with a growth of the last term in 
each equation. As discussed in Sec. IV, this is precisely 



defined trajectòries out of an ensemble spreading out in 
phase space. 

At the level of these second cumulants, the exact form 
of the damping is, thus, immaterial. The fact, however, 
that we derived these from the stochastic master equa- 
tion (|J) gives us not only a theoretical 'unraveling' of 
the master equation, but also provides a physical mean- 
ing to each trajectory. Furthermore, it guarantees that 
the underlying evolutions of density functions are com- 
pletely positive, and, therefore, not only does the co- 
variance matrix stay positive, but also these equations 
can be completed into a hierarchy of cumulant equa- 
tions that automatically satisfy the appropriate reality 
conditions. Furthermore, as discussed in Sec. [ÏÏ|, such 
a measurement process leads to the unavoidable noise 
proportional to Ti 2 k apparent on the right hand side of 
Eq. (||); in contrast, in the classical discussion in Sec. |í"v| , 
the corresponding noise term is unrelated to the mea- 
surement process and can even be set to zero. In fact, 
the truncation to the second cumulants implies that no 
truly quantum effects of dynamics come into play p5[ | in 
our approximation, and the quantum scale h appears in 
our equations purely from this information-disturbanec 
consideration pq . 

We will make two self consistent approximations in 
order to examine in what regime classical dynamics 
emerges. The first is to truncate the power series in 
F to second order. The second is to neglect third and 
higher cumulants in the equations for the second cumu- 
lants. An examination of the equations of motion for the 
third cumulants appearing in Eqs. (||) to ( |ï(ï| ) shows that 
indeed these are damped by the measurement, again with 
damping coefficients proportional to k. The fact that the 
wavefunction stays close to Gaussian is also borne out by 
numerical simulations pij . 

Setting the third cumulants to zero in the equations 
for the second cumulants, we solve for the stable steady- 
state: 



V ss = 



2C 



Xp 



mk 

mV* s (kC x s p -d x F). 



(15) 
(16) 



CS; 



d x F , s fd x F 

— r~ +sgn(m)W — r - 



h 2 

4:1] 



T-- (17) 



the way in which classical measurement also seleets well 



where sgn(m) is the sign of m which we shall henceforth 
take to be positive [^7j and d x F is taken to be evaluated 
at a typical point in phase space. Now, there are three 
conditions that must be satisfied in order for the clas- 
sical limit to be obtained. First, localization such that 
V x d x F(x) <C 2F(x), as discussed above, must be main- 
tained, second, the noise introduced by the measurement 
should be negligible compared to the classical motion, 
and third, that the measurement record should follow 
the motion of the position with sufficicnt aceuracy. 

Before examining these conditions in turn, two points 
are in order. First, localization and low noise are not 
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really independent constraints: in fact to provide effec- 
tive damping for the covariance matrix, the noise has to 
increase with increasing width of the state. Conversely, 
noise also cffects a spread in phase space of any uncertain 
state, especially near unstable points. It is convenient, 
however, to treat the direct effect of the finite width of 
the state on the deterministic evolution in a nonlinear 
potential as a question of localization, and the rest as a 
question of low noise on the trajectòries. 

As a second point, it is important to emphasize that 
the deviations of the quantum trajectòries from the clas- 
sical ones can be diffcrcnt in different parts of the phase 
space. Nevertheless, in most experimental situations, 
'classical' quantities evaluated are of similar orders of 
magnitude almost everywhere, and, so, we shall ignore 
these differences and consider them evaluated at a 'typi- 
cal' point around the trajectory in question. 



A. Localization 

We start by noting that for the deterministic part of 
the equations of motion for the quantum mean vàlues x 
and p to match the classical equations of motion, we need 



(F(X)) = F(x) + \v x d 2 x F{x) + 

to very closely approximate F(x). That is, we need 
_ d 2 x FV x 



2F 



«1. 



(18) 



(19) 



Replacing V x with its typical steady state value 
[Eqs. ( |15| , [l7| )1, we see that r is the positive solution to 



2mk 2 F 2 



■r 2 -(d x F) = J(d x F) 2 + 



2 J.2 



h z k 

if] 



(20) 



Since this equation implies that r is a monotonically de- 
creasing function of fe, Eq. ( |Ï9| ) can provide a lower limit 
for k. We examine this possibility in the following dis- 
cussion. 

Due to the positivity of its right hand side, Eq. j20| ) 
implies that at the unstable point d x F > 0, we must have 

2mk 2 F 2 r 2 > (d 2 x F) 2 d x F . (21) 

This alone means that to have r <C 1, it is necessary that 



2 „ {d 2 x F) 2 \d x F\ 



iF 2 



(22) 



Squaring Eq. ( p0| ) one sees that r is the algebraically 
largest solution of 



F 2 r 4 \ IQrjmF 2 



(23) 



For this solution to be small would generically require 

k > ; Z_ , 24 
except in the typical case of small nonlinearity 

n 

when the width stays small at the stable points indepen- 
dent of the value of k. 

Since, when the nonlinearity is large enough to violate 
Eq. ( p5[ ) , Eq. (p4j) is stronger than Eq. (p3) , we can sum- 
marize these results as follows: If the nonlinearity, char- 
acterized by d x F, is sufhciently weak to satisfy Eq. (EH), 
then, at the unstable points (d x F > 0) one needs 



8rjk > 



(d 2 F) 2 \d x F\ 



2mF 2 

In the case of strong nonlinearity, we need 

(d 2 x F) 2 h 



8rjk > 



A^/rjmF 2 ' 



(26) 



(27) 



to hold at all points. 



B. Low Noise 

In this subsection we consider the noise component of 
these equations. In the classical limit the effect of this 
noise must be negligible on the scale of the determinis- 
tic dynamics. To compare the random noise with the 
deterministic dynamics, we need to average over an ap- 
propriate time scale: during a time T upon which the 
dynamics is effectively linear, the noise dW provides an 
rms contribution of \ZT. We will define 'low noise' to 
mean that the noise contribution on this time scale is 
small compared to the deterministic contribution. The 
time scales upon which the dynamics is linear for x and p 
are those in which the terms in the respective equations 
do not change appreciably, and we will use the deter- 
ministic terms to obtain these time scales. The deter- 
ministic motion for x is driven by p/m, so appreciable 
changes occur when the change in momentum, Ap, is of 
the order of p. The resulting time scale for this change 
is T x ~ \p/F\. The deterministic motion for p is driven 
by (F(X)) w F. Changes in F are due to changes in 
x: in particular AF w d x FAx. Hence the time scale 
for changes in F is T p ~ (m\F\) / (\pd x F\). Demanding 
that the change in x and p from the noise is small com- 
pared to that due to the deterministic motion in these 
time intervals gives the two inequalitics 



If 3 



vH, « - v ~ , rl • 

m 1/ Tïl\r 



\/~kCxp 



< FJT V = 



\F*\ 



\pd x F\ 



(28) 
(29) 
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We will now examine these two inequalities in turn. 

Considering the first inequality, and replacing V x with 
its typical steady state value, given above, we have 



(30) 



where E = p 2 /(2m) is the typical energy of the system. 
Noting that Hlpl/lí"! has units of action, to simplify the 
following analysis we will now define a dimensionless ac- 
tion s by sfr = E\p\/ (4|.F|). Conceptually s may be iden- 
tified with the typical action of the system in units of h. 
Using Eq. (|ï(]) to write C x s p in terms of k the inequality 
becomes 



'd x F\ h 2 „ d x F 
k J 4/7 k 



(31) 



The positivity of the left hand side immediately gives us 
the condition 



d x F 

hk 



< As. 



(32) 



Now squaring both sides of Eq. ( |3Ï| ) , and rearranging, we 
obtain 



< 167Ts 



2„2 



8hsd x F 



4?7 k 
Because of Eq. ( |32] ) , this condition reduces to 

1 

s > 



8VV 



(33) 



(34) 



cxccpt at the unstable points (d x F > 0) where one re- 
quires, in addition, 



hk 32ns 

> 



1 



d x F 64r)s 2 - 1 2s 



(35) 



where the approximate equality is implied by the inequal- 
ity in Eq. Q. 

We now consider the inequality given by Eq. ( |29| ) . Re- 
placing C xp with its typical steady state value, and per- 
forming some rearrangements we obtain 




£ « 4r? Is' - sgn{d x F) 
where for compactness we have written 

hk 



(36) 



hs' 



\d x F\ 

mF 2 \F\ 

(d x F) 2 \p\ 



(37) 



Here s' is a dimensionless quantity, which we will once 
again take to be an estimate of the typical action of the 
system in units of h. For d x F > 0, the condition 



f < 4r? Is' - 




(38) 



is satisfied whenever 



16 , „ , 

— < Ç < 2-qs , 



(39) 



whereas, for d x F < (the condition is not useful when 
d x F = 0), it is sufficient that 



hk 

Wf\ 



(40) 



Collecting all the inequalities in this subsection, (i.e. 
Eqs. @, (g), @, © and ©), we find that they 
are all implied by 



r)s 4 



(41) 



where s = min(s, s'). 



C. Faithful Tracking 

In the previous subsections we have been considering 
the motion of the centroid of the quantum wave packet, 
(x,p). This centroid represents the observer's true best- 
estimate of the mean value of position and momentum 
at the current time, given the measurement record. To 
obtain this best estimate the observer must know the 
dynamics of the system, given by Hem, and then integrate 
the full stochastic master equation, where the correct dW 
is obtained continuously from the measurement record. 

In practice, it is often merely the measured value of 
position which is taken as the estimated value. Hence, 
we need to find conditions under which this value tracks 
the true best-estimate with sufficient aceuracy. Since the 
measurement record in our formulation contains white 
noise, the simplest way to model a realistic macroscopic 
measuring apparatus is to low-pass filter, or band limit 
the measurement record to obtain the continuous esti- 
mate of the position (this is equivalent to making the rea- 
sonable assumption that all real measuring devices have 
a finitc response time). This is achieved by averaging the 
measurement record over some finite time At. To obtain 
an aceurate estimate, At must be short compared to the 
dynamical time scale of the system. 

If we assume that the change in x over time At is 
negligible, then the error in the estimate of x resulting 
from averaging the measurement record y(t) over Aí is 



a T (x) = {SrikAt)- 1 / 2 . 



(42) 



Hence, if to accurately track a classical dynamical system 
we require a spatial resolution of Ax, and a temporal 
resolution of Aí, then we must have 



8r]k > 



1 



At(Ax) 2 



(43) 



We also note, however, that in the observation of classi- 
cal systems, classical estimation theory is, in fact, often 
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used to obtain the classical equivalent of the quantum 
best-estimates provided by the SME (Eq. (||)). Such a 
procedure is most often used in classical feedback-control 
applications. In the classical limit, therefore, such classi- 
cal estimation procedures must work effectively, and we 
will verify this in the next section. 



D. Summary 

We have now derived a set of inequalities which, whcn 
satisfied, lead to the emergence of classical mechanics. 
Consider first the inequalities which come from the lo- 
calization condition. In the macroscopic regime, which 
applics to common mechanical devices one would build 
in the laboratory, the right hand side of inequality (|25| ) 
is extremely large compared to the typical nonlinearity. 
Consequently this inequality is satisfied, and the rcsult- 
ing condition for k is given by (p6|). Note that h does 
not appear in this inequality. In fact, this is actually a 
classical inequality, similarly required for classical con- 
tinuous measurement on classical systems. In that case, 
the observer's state of knowledge of the system is given 
by a classical probability density in phase space, and this 
evolves as the system evolves and as information is con- 
tinuously obtained. 

If the system is sumciently small, and the nonlinearity 
sumciently large on the quantum scale so that inequality 
( p5| ) is not satisfied, then the condition for k is replaced 
by inequality (0). This does contain h, and is, therefore, 
a uniquely quantum condition. It appears due to the 
unavoidable quantum noise which affects the dynamics 
strongly if the nonlinearity is large on the quantum scale. 

The left inequality in Eq. (JO]), again is a classical con- 
dition (the h arises because we chose to measure the ac- 
tion s in units of fi) : it reflects the observation that if the 
measurement does not localize the motion, the state es- 
timate changes from moment to moment essentially ran- 
domly, or in other words, the noise is large. The right 
hand inequality in Eq. ([4l]) is the direct effect of the 
irreducible noise coming from the measurement process 
and is thus a quantum effect. Together, as the action 
increases, these low noise conditions put ever decreasing 
constraints on the required measurement strength. 

The faithful tracking condition is once again purely 
classical, in that it also applies to classical observation. 
It is simply the condition on the accuracy of the measure- 
ment so that the measurement record itself, as opposed to 
the estimated state, accurately tracks the motion of the 
system from which the localization condition is derived. 

It is worth noting that the above inequalities also de- 
tcrmine the regime in which múltiple observers agree on 
the motion of an object, which is clearly an important 
property of the classical limit. As discussed in Section ||, 
múltiple observers can be taken into account by giving 
each observer, i, a value of r\ — rji such that J^. rji < 1, 
and giving each a different noise realization, dWi. Fur- 
thermore, it is clear from the derivation of the stochastic 



master equation 1 23 , g4| that the state conditioned by the 
measurements made by all of the observers is narrower 
than and consistent (in probability) with the state esti- 
mate of each observer; ipso facto, the estimates of the 
different observers must agree within errors. Since the 
conditions derived in this section can be satisfied with 
rj < I (even with ïj< 1), and since these imply localiza- 
tion and accurate tracking of the measurement record, 
under these conditions all observes will agree upon the 
motion of the system to errors small on the classical scale. 



IV. CLASSICAL ESTIMATION IN THE 
CLASSICAL LIMIT 

When a classical system is subject to noise and contin- 
uous observation, a classical theory of continuous state- 
estimation may be developed to describe the continu- 
ous acquisition of information regarding the system p8| . 
Consider an observed classical system whose dynamics is 
given by 



dx 
dp 



p/m 



dt - 



F c (x 

with measurement record 

dy c = xdt - 







dV 



(44) 



(45) 



i.e., we consider a system with purely additive momen- 
tum noise being observed continuously and with random 
errors. Here, dW p and dV are Wiener noises with dV 
possibly correlated with the dW p , and g p and g m are 
positive real numbers. Then the evolution of the state 
of knowledge of the observer, described by a probability 
density P{x,p,t) obtained by averaging over dW p and 
conditioning by dy c , is pq , p9| 

dP = [-(p/m)d x -(F c - g p d p )d p ]Pdt 

+Vg^{x - (x))PdW , (46) 

where dW — ^Jg m (x — (x))dt + dV, and turns out to be 
a Wiener noise, uncorrelated with the conditional prob- 
ability P. Note that we can then write the measurement 
record for the classical measurement as 



dy c = (x}dt + 



dW 
\[9^ 



(47) 



and we see that this can be viewed as directly analogous 
to the quantum measurement record. The equations of 
motion for the classical best estimates (x) c and (p) c , and 
the second order moments are 



d(x) c = ^dt 
m 



WxdW , (48) 
d(p) c = (F c (X))dt + ^C xp dW, (49) 
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and 



dV x 

dVr, 



m 



9mVx 



dt 



dC 



[2g p - g rn C 2 xp + 2d x FC xp ] dt 
+d x FK xxp dt + y/g~^K xpp dW , 
1 




x p 



Vp — g m V x C xp + d x FV x 



dt 



+ -d x FK xxx dt + ^/g^K xpp dW . 



(52) 



Identifying g rn = k and g p = fi 2 k, we see that these equa- 
tions are identical to the quantum equations governing 
the continuously estimated state [Eqs. (P|)-(|ÏC|)], and the 
only way that quantum mechanics enters is in enforcing 



fí 9m < 



Even though this is the case, it should be 



noted that when the potential is nonlinear, the equations 
of motion for the third and higher cumulants are not the 
same in the quantum and classical cases, so in general the 
evolutions of the classical and quantum estimates differ. 
In the classical limit, however, the conditional probabil- 
ity or the state in the quantum case, is Gaussian to a very 
good approximation so that the third cumulants can be 
set to zero, and as a result they no longer feed into the 
equations for the second order cumulants. Consequently, 
the evolution of the classical best estimates and second 
cumulants are identical to the quantum estimates for the 
same measurement record, and as a result classical esti- 
mation may be used to track dynamical systems in the 
classical limit. 



V. NUMERICAL EXAMPLES 

In this section we provide numerical support for the 
arguments in the previous section. We present two exam- 
ples, and show that under the conditions derived in the 
previous sections, the quantum wave packet remains lo- 
calized, the evolution of the centroid follows the classical 
motion with negligible noise, and both the measurement 
record (suitably band limited) and the classical state- 
estimate accurately track the motion of the system for 
each of a set of observers. 

To derive the equation of motion for the wavefunc- 
tion of the continuously observed system, assuming 
N observers, one can first write down the Stochastic 
Schròdinger equation for the unormalised wavefunction 
for a single observer, making N measurements. If the 
interaction strength for measurement i is r^fc, then this 
is 



\{iH(t) + HkX 2 )dt 



N 



^rjikdn 



i=l 



(53) 



where the record for each measurement is given by 

dW,. 



dr % = (X)dt + 



y/8r^k 



(54) 



Now we let each observer have access to just one of the 
measurement records. In addition, we choose J2i Vi — lj 
so that r\i represents the fraction of the total measure- 
ment interaction strength fe used by each observer. The 
evolution of the state-of-knowledge for any particular ob- 
server (who only has access to her measurement record) 
can be calculated by averaging over the noise realiza- 
tions for all the other observers while keeping the mea- 
surement record for the observer in question fixed. The 
resulting equation of motion for the state-of-knowledge 
of observer i, given the measurement record dri gener- 
ated by the stochastic Schròdinger equation (^3|) is the 
stochastic master equation |23|, H 



dp = — — [H,p]dt 
n 



k[X, [X,p}}dt 



-([X,p} + -2pTr[pX])y/2 V ~kdV l (55) 



where 



dVi = y/8rük(dn - Tr[pX]dt). 



(56) 



Notc that this is, in fact, just Eq.(^|), because as far as 
observer i is concerned, all the other observers are sim- 
ply gathering part of the environment to which i has no 
access. In addition, the fractions r\i are the respective 
measurement efhciencies. 

To simulate múltiple observations on a given sys- 
tem, we first integrate the stochastic Schròdinger equa- 
tion which generates a set of measurement records, one 
for each observer. We then integrate the correspond- 
ing stochastic master equations using the measurement 
record for each observer. The state-of-knowledge of each 
observer over time can then be compared to the 'actual' 
evolution of the system state vector given by Eq. (|53|). 
The stochastic Schròdinger and master equations were 
integrated in time using a spectral split-operator method. 
Since the classical limit is obtained when the extent of 
the wavefunction is small compared to the range of mo- 
tion of the centroid, the algorithm is designed so that 
the computational grid follows the wavefunction in both 
position and momentum space, and this is crucial for ef- 
ficient computation. 

In treating systems of different sizes and actions it is 
convenient to choose units for the system variables to 
keep the numerical value of the action close to unity. 
Due to this system dependent choice of units, the fixed 
quantity h has a system dependent numerical value; and 
indeed we expect the classical limit when fi <C 1 in these 
units. This is what we demonstrate below. 



9 




t 

FIG. 1: The Standard deviations of the state estimates for 
each of the three observers for the Duffing oscillator, plotted 
over a duration of t = 5. Sòlid line: observer with rj — 0.5; 
dashed line: observer with 77 = 0.3; dash-dot line: observer 
with 77 = 0.2. 



A. The Duffing oscillator 
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t 

FIG. 2: Plot of the error Standard deviation (blue) and the 
difference between the estimated and true means for a sin- 
gle noise realization (red) for the simulation of the Duffing 
oscillator using the stochastic master equation with r\ — 0.5. 
The green and purple curves plot the same quantities for the 
Gaussian estimator. 



The Duffing oscillator is a sinusoidally driven double- 
well potential, with Hamiltonian 

H{t) = — + BX 4 - AX 2 +AX cosí Lüt) . (57) 

Wechoosem = 1,A=1Q,B = 0.5, A = 10 and ui = 6.07. 
At times when the driving is zero, this puts the minima of 
the two potential wells at ~ ±3.2, with a central barrier 
height of 50. We choose h = 10~ 5 and k = 10 5 , which is 
sufhcient to satisfy the inequalities derived in Section III 
for all but tiny vàlues of 77, and therefore puts the system 
in the classical regime. We now evolve the system with 
three observers, and set their measurement eficiències to 
be 0.5, 0.3, 0.2 respectively. We first calculate the posi- 
tion variance Vx of the wave-function given by evolving 
Eq. (p3|)), and verify that this remains sufficiently small. 
Running the simulation for a duration of t = 5, the max- 
imum value of y/Vx is 2.7x 10 -3 , and the rms value over 
the evolution is 1.4 x 10~ 3 . The localization condition is 
therefore well satisfied, and an inspection of the evolution 
of the centroid shows that the noise is indeed negligible. 
Showing that the evolution is indeed the classical evolu- 
tion is more nontrivial since the system is chaotic: any 
small difference in the noise on two trajectòries will cause 
them to diverge rapidly, and one cannot therefore sim- 
ply compare the trajectory to the equivalent noise-free 
classical trajectory. In Ref. fïij , the classical dynamics 
was verified by comparing the stroboscopic map and the 
largest Lyapunov exponent obtained from the quantum 
evolution and their classical equivalents. Here we calcu- 
late the continuously estimated state, both quantum and 
classical, for the different observers, and show that these 



agree, and agree between observers. 

We now calculate the quantum state-estimate for each 
observer, obtained by integrating Eq. J55|), and compare 
this with the classical (Gaussian) estimate for each ob- 
server, obtained by integrating Eqs. (fÏ8|)-([5"2|). In Fig- 
uro Q we plot the uncertainty in position (characterized 
by y/Vx = y/Tr[pX 2 } - Tr[pX] 2 ) for the quantum state 
estimated by each observer over the duration of the run. 
All these remain small. The rms of yT4 for each ob- 
server over the duration of the run is given in Table Q. 
The evolution of the uncertainty in position for the Gaus- 
sian state-estimate is essentially identical to the quan- 
tum estimate, and the rms of \fVx for this estimate is 
also given in Table |. Note that the position variances 
for each observer are, as expected, larger than the vari- 
ance of the wave-function calculated using the stochas- 
tic Schròdinger equation. Since the solution to Eq. ( p33] ) 
can be viewed as an 'unraveling' of the stochastic master 
equation in Eq. (|55|), the difference between the variances 
of the 'true' state estimate from the former and the indi- 
vidual observers' state estimate from the latter provides 
an estimate of the amount by which their means differ, 
averaged over noise realizations for all the observers. We 
will refer to this, for want of a better name, as the error 
variance, and its square root as the "error Standard devi- 
ation" . The rms value of this error Standard deviation is 
1.2 x 10~ 3 , 1.7 x 10~ 3 and 2.1 x 10~ 3 for observers with 
i] = 0.5, 0.3 and 0.2 respectively for both the stochastic 
master equation and Gaussian simulations. In Figure ^ 
we plot the evolution of the error Standard deviation, and 
also the actual difference between the estimated mean 
and the mean of the wave-function, for both the master 
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TABLE I: rms Standard deviation of state-estimates for the 
Dufhng oscillator and the rms deviation of the averaged mea- 
surement record. 



Observer's 77 


77 = 0.5 


77 = 0.3 


77 = 0.2 


Quantum 
Gaussian 
Averaged Record 


1.9 x 10" 3 
1.9 x 10 -3 
8.2 x 10" 3 


2.3 x 10" 3 
2.3 x 10" 3 
9.3 x 10" 3 


2.6 x 10" 3 
2.6 x 10" 3 
1.1 x 10" 2 



equation simulation and the Gaussian estimator. The 
small difference between these last two is most probably 
due solely to the fact that the mean of the computation- 
ally intensive master equation simulation has not com- 
pletely converged at the value of the time step employed. 
(This difference is not seen for the computationally sim- 
pler delta-kicked rotor system discussed in the next sub- 
section.) 

Each observer may also track the position simply 
by averaging her measurement record over a suitable 
time period (i.e., by low pass filtering the measurement 
record). Naturally this period should be as long as pos- 
sible so as to filter out the noise, but short enough so as 
not to filter out the deterministic motion. For this sys- 
tem we use a time period of 2.5 x 10~ 2 for the filtering. 
The average rms deviation of this estimate from the mean 
position of the wave-function is also given in Tablc | for 
each observer. From this we see that all observers can 
effectively track the motion of the particle (up to an er- 
ror in position of about 10~ 2 ) using their measurement 
records directly. 



B. The delta-kicked rotor 

The delta-kicked rotor obeys the Hamiltonian 

p2 °° 

HOt) = — — T K cos(X) V SOt - n) . (58) 
2m z — ' 

n=0 

It is, thus, a free particle, which experiences regular kicks 
from the potential of a nonlinear pendulum. For a wide 
range of parameters, the quantum behavior of this system 
(by which we mean the evolution of the closed system) 
is very different from the classical motion. In particular, 
after a few kicks the average energy of the closed clas- 
sical system increases linear ly with time. In the closed 
quantum system, however, the average energy reaches a 
maximum value and after that point remains fairly con- 
stant. This is termed dynamical localization. We now 
simulate the evolution of the observed wave-function for 
this system, with the same vàlues of fi and k as we used 
for the Duffing oscillator, and with the same three ob- 
servers. For the system parameters we will choose K = 10 
and m = 1, and integrate for a time period of 30 kicks. 
First we check the localization of the wave-function given 
by integrating the stochastic Schròdinger equation, and 
find that the average value of \JVx is 2.1 x 10~ 3 , and the 
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FIG. 3: The average kinetic energy for the delta-kicked rotor 
as a function of time. The classical value is obtained by aver- 
aging over 10,000 trajectòries. The observed quantum value 
was obtained by averaging over 1000 trajectòries. 
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FIG. 4: Plot of the error Standard deviation (green) and the 
difference between estimated and true means for a single noise 
realization (red) for a simulation of the delta-kicked rotor us- 
ing a stochastic master equation with 77 = 0.5. The results for 
the Gaussian estimator are indistinguishable on this seale. 



maximum value obtained during the run is 3.2 x 10 . 
We check that the mean energy is indeed behaving in a 
classical fashion by averaging this energy over many re- 
alizations, and comparing this to the classical value. In 
Figure [|we plot the average energy of the observed quan- 
tum system, using h = 0.1 and k = 10, along with both 
the classical result and the quantum result for h = 0.1. 

We next compare the position uncertainties in the 
state-estimates of the different observers, as above for 
the Duffing oscillator, and present these results in Ta- 
ble pi The Gaussian estimator agrees with the stochastic 
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TABLE II: Average deviation of state-estimates for the delta- 
kicked rotor and the rms deviation of the averaged measure- 
ment record. 



Observer's 77 


■q = 0.5 


77 = 0.3 


r) = 0.2 


Quantum 
Classical 

Averaged Record 


2.9 x 10" 3 
2.9 x 10~ 3 
8.6 x 10" 3 


3.6 x 10" 3 
3.6 x 10" 3 
1.0 x 10" 2 


4.3 x 10" 3 
4.3 x 10 -3 
1.1 x 10" 2 



master equation, and the uncertainties are small, so that 
the observers effectively all agree on the motion. The 
averaged measurement record also tracks the motion ef- 
fectively. The rms valuc of the error Standard deviation 
is 1.9 x 1CT 3 , 2.9 x 1(T 3 and 3.7 x 1CT 3 for observers with 
?; = 0.5, 0.3, and 0.2 respectively. In figure |^ we plot the 
evolution of the error Standard deviation for the observer 
with 77 = 0.5 and also the actual difference between the 
estimated mean and the mean of the wave-function for a 
single realization of the stochastic master equation sim- 
ulation. The equivalent plots for the Gaussian estimator 
are virtually indistinguishable. 

To conclude, we see from the above simulations that (i) 
in the classical regime the full quantum state-estimation 
reduces to Gaussian state-estimation, and hence classi- 
cal state-estimation may be used, (ii) even without the 
use of true (and therefore optimal) state-estimation, low 
pass filtering of the measurement record alone provides 
adequate tracking of the system, and (iii) since the er- 
rors in the respective estimates are small, all observers 
effectively agree upon the motion of the system. 

VI. CONCLUSION 

The emergence of classical dynamics remains a central 
issue in understanding the predictions of quantum me- 



chanics, especially now that experiments are becoming 
available to probe this transition directly jï(|. In this 
paper, by dcriving general incqualities which determino 
when classical mechanics will emerge, and by providing 
numerical examples, we have presented very substantial 
evidence that quantum measurement theory provides a 
completely satisfactory answer to the question of how 
classical mechanics, and hence classical chaos, emerges 
in a quantum world. In doing so we have shown in some 
detail how the mechanism for this transition can be un- 
derstood as a result of localization and noise suppression 
in the classical regime. 



While the emergence of classical dynamics for a single 
motional degree of freedom now appears to be well un- 
derstood, the quantum to classical transition as yet holds 
many unanswered qüestions. What happens, for exam- 
ple, to the dynamics of a system as it passes "through" 
the transition? How do systems behave when they are 
neither fully quantum nor fully classical? For example, it 
is known that the delta kicked rotor demonstrates a com- 
plex behavior in the transition region p0| . Further qües- 
tions include how classical dynamics emerges for other 
degrees of freedom, such as spin, and what happens, for 
example, when spin and motional degrees of freedom are 
coupled? Must all the subsystems have a large action (we 
note that this has recently been investigated, see pit), 
and must all the degrees of freedom be continuously mea- 
sured, or will a subset sufhce? For a spin system, must 
one measure all the components of spin, or will a single 
component suffice? Fortunately we are now at the point 
where one can not only pose these qüestions, but expect 
that sòlid answers will soon be forthcoming. 
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p — > p* is an invariance of Eq. ^| Thus our results also 
apply to the negative mass case mutatis mutandis. 
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